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The focus of this paper will be on the development of a photovoltaic (PV) 
pumping system that employs a DC-DC boost converter as an intermediate 
power conditioning unit in a centrifugal water pump powered by a brushless 
DC (BLDC) motor supplied by a solar photovoltaic panel (SPV). This work 
describes a particle swarm optimization (PSO) maximum power point 


tracking (MPPT) method compared with the gray wolf optimizer (GWO) 
method that are closed-loop with a sliding mode controller (SMC) in order to 
control the operational point of a solar PV in partial shading conditions (PSC), 
the DC-DC boost converter provides a low starting current for the BLDC 
motor, and the voltage source inverter (VSI) is used to perform electronic 
commutation of the BLDC motor. Furthermore, the BLDC motor's speed is 
regulated by the changeable DC link voltage. The system under consideration, 
as well as the control approaches employed, were created in the 
MATLAB/Simulink environment. 
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1. INTRODUCTION 

Rural and remote areas where there is no access to electricity networks, suffer from an increased lack 
of water. This shortage has led to a growing interest in the use of photovoltaic (PV) generators as an energy 
source. Indeed, the practical and cost-effective solution to the problem of water scarcity is to implement a 
stand-alone PV pumping system. Motors absorb more than 40% of the energy consumed in driving systems. 
As aresult, motors are critical in applications requiring electric drives, such as pumping water in a solar system 
[1]. Induction motors AC first drew a lot of attention for water pumping systems [2], [3] due to its simplicity, 
low construction cost, and low maintenance. However, as speed increases, brush friction necessitates more 
precise control, reducing useful torque [4]. The hunt for a superior substitute for induction motors led to the 
creation of brushless DC (BLDC) motors. In this regard BLDC motors have become popular in the automotive, 
robotics, and medical industries in this respect. These motors have several advantages over other types of 
motors, including silent operation, extended life, excellent efficiency owing to smooth speed control, enhanced 
power density, and reduced electromagnetic interference [1], [5]. 

The power efficiency of the gross premium valuation (GPV) has a major influence on the performance 
of a PV pumping system. Because PV systems are non-linear, the use of an maximum power point tracking 
(MPPT) control rule is required to assure maximum power transfer to the load [4]. Several MPPT approaches 
are available in the literature to reach this maximum power point (MPP) [6]-[8]. Among the most commonly 
used traditional approaches are the perturb and observe (P&O) method, the incremental conductance (INC) 
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method, the fractional short circuit, and the fractional open voltage circuit [7], [8]. However, because the MPP 
varies with sunshine and season, it is difficult to keep this sort of MPP tracker running. There was a 
considerable loss of electricity. Due to the difficulty in detecting the precise MPP when the P-V characteristics 
include several peaks [9], [10], the system has become less efficient under changing irradiation conditions [11]. 
To solve this problem, many improved MPPT algorithms, fuzzy logic, and artificial intelligence-based systems 
have been created, in order to assess the system's dynamic reactivity, and provide better MPP monitoring. This 
work describes a particle swarm optimization (PSO) method compared with the gray wolf optimizer (GWO) 
method that have closed-loop with a sliding mode control (SMC) for a weather-aware PV water pumping 
system. SMC is recommended as a particularly suitable approach to cope with variations in solar radiation for 
PV system operation [11], [12]. In general, the selection of DC-DC converters is critical in defining the best 
performance of BLDC drives, and as a result, enhancing the output of the renewable energy system [3]. In the 
PV system, various types of DC-DC converters, such as buck, boost, buck boost, cuk, and single ended primary 
inductor converter (SEPIC) [3], [8], [9], [13], [14], are implemented with the following considerations in mind: 
high output voltage gain, low input current, low current ripples, higher efficiency, minimal cost, higher 
efficiency, high maximum power flow, and easy system drive [15]. Figure 1 shows a schematic design of the 
investigated water pumping system. This system is made up of PV solar panels, a boost DC-DC converter, a 
three-phase inverter voltage source inverter (VSI) operated BLDC motor and a centrifugal pump. 

The following is an excerpt from this article: in section 1, there is an introduction. Section 2 presents 
the studied system's configuration (solar photovoltaic panel (SPV) network, DC-DC boost converter, BLDC 
motor, and centrifugal pump). Section 3 presents the studied system's control (PSO and GWO algorithms and 
SMC). Section 4 presents the results and analysis using the MATLAB/Simulink environment and conclusion 
in section 5. 


PV Array 


BLDC Motor 
hg} bdh, 
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6 Switching Centrifugal Pump 


Sequence 
Generator 


Figure 1. Solar pumping system design 


2. SYSTEM DISCRIPTION 
2.1. Photovoltaic panel 

PV module efficiency varies with the quantity of light received. Figure 2 depicts the corresponding 
circuit model of a solar cell [9] and is characterised by (1)-(4): 


lp = lpn — Is jee) 5 | pan (1) 
lpn = [pno + Ju(T - Teol È 2) 
Is = lso (Z) e ae ae?) 3) 
Iso = oo (4) 


eAKNsTsc—1 


With Ipy and Vp, are respectively the PV panel current and voltage, N, is the number of PV module cells, K is 
the Boltzmann constant (1.3806503e—23 J/K), T is the cell temperature in Kelvin, q is the electron charge 
(1.60217646e—19 C), Iso is the saturation current at Tse , Inno is the photo-current measured under the standard 
conditions E,,, Tse (Esc =1,000 W/m? and T,, =25 °C), and J; is the temperature coefficient of the short circuit- 
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current. The key parameters for utilizing this model in simulation are shown in Table 1. To illustrate the 
behavior of the PV model, the current-voltage and power-voltage curves of the model for different irradiation 
values at constant temperature (25°) are shown in Figure 3. 

PV module efficiency varies with the quantity of light received as presented in Figure 4(a). As a result, 
it is critical to investigate the qualities of a PV module under various irradiation scenarios, such as uniform 
irradiation and partial shading conditions (PSC), where multiple local and global MPPs may be seen on the P-V 
characteristics curve [6], [10] as seen in Figure 4(b). The use of bypass diodes linked in parallel to each PV module 
during PSCs reduces the possibility of a hot spot, in which the shaded module works as a load rather than 
generating power [16]. The top peak in Figure 4(b) for the PV array is the global peak, which is successively 
equivalent to 188.53 W for profile 1 and 131.89 W, while the others are local peaks. 


Table 1. Parameters of the PV array Tio 

Parameters Nomenclature Values S 5 3 
Voltage at open cct Voc 21.7 V rs] e A P E a aa ao ag 
Short cct-current Ise 7.3 A Voltage (V) 
Max power voltage Vinpp 32.82 V 100 y ' 
Max power current Lapp 8.07 A 
Photo-current Iph 8,6307 A 
Saturation current I, 1.4176e-%°A 2 á 6 i 10 2 na 
Diode ideality A 0.99132 ee 
Series resistance Rs 0.098625 Q . 
Shüñtresistáncé Ren 82.1161 Q Figure 3. The current-voltage and power-voltage curves 
Bandgap’s energy Egap 1.12eV of the model 
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Figure 4. Under PSC: (a) structure of the PV panel and (b) the I-V and P-V characteristics of the PV panel 


2.2. DC-DC boost converter design 

Choosing DC-DC converters in general is the best approach for determining the optimal performance 
of BLDC motors and therefore optimizing solar system production. Indeed, the SMC is designed to follow the 
voltage reference V,¢¢ by adjusting the duty cycle of this DC-DC boost converter. Figure 5 is an example of a 


boost converter, and its average model is given by (5): 
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dVpv 1 
T = = (pv ~ I) 
TL = E (Vy, — (1 — Vier) Ve) (5) 
dt  L\'PY réf/ Ys 
dVs 


1 
F zala = Ver) = Is) 


The relation between the input and output voltages is shown in (6): 


a (6) 


Figure 5. The boost converter circuit 


2.3. Modeling of brushless DC motors, commutation, and voltage source inverter switching 

BLDC motor is also called electronically commutated motor [7], [8], [13]. This electric actuator is said 
to be the best for solar pumping in Figure 6. The main purpose of electronic commutation is to operate a VSI in 
such a way that the DC current drawn from the DC link capacitor is symmetrically placed over 120° in the center 
of each phase. To achieve this, the rotor position must be monitored at each switching time, i.e., every 60 electrical 
degrees or six times in an electrical cycle. To do this, the position information is taken by a position sensor, often 
a hall sensor [17]. When the magnetic poles of the rotor approach the hall sensors, they create a separate set of 
three hall signals (hy, hz, and h3). These signals are digital in nature and are logically transformed into six 
switching pulses (S1, S2, S3, S4, S5, and S¢) that are used to control the VSI switches [14], [18]. 
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Figure 6. Configuration of the proposed SPV array and boost converter powered by a BLDC motor for a 
water pumping system 


The comparison of the duty cycle generated by the MPPT/SMC assembly with a high frequency 
sawtooth carrier wave (to keep the symmetry) gives us a high frequency pulse width modulation (PWM) pulse. 
A logic AND gate is used for fundamental frequency modulation [8], [14]. As an example, the conduction 
states of two switches (s4 and s4) of a VSI supplying a BLDC motor are shown in Figure 7. 


Figure 7. A BLDC motor powered by a VSI during the conduction states of (sand s4) 
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The MPPT algorithm generates a duty cycle for the VSI, which controls the link voltage Vpc. Then, as 
the atmospheric conditions change, so does the power supply by the PV and the speed of the BLDC 
motor [13], [18]. The role of power conditioning is to maximize the energy transfer between the PV array and the 
pumping unit. The power conditioning is provided by a DC/AC inverter. The three-phase inverter voltages 
(Va Vp, and V,) are represented by (7)-(9) [15]. It should be noted that the BLDC motor stator phase currents (ig, 
ip, and i,), flux links (Ag, Ap, and A,), and back emfs (ea, €p, and ec) have trapezoidal waveforms [13], [19]. 


Va = Raia + Ag+eq with Ag = (La — Lm)” 
Vp Bs Rpip + Àp + ep with Ap = (Lp a Lin ae (7) 


Vo = Reic + Ate, with — Ae = (be — bm) <2 


Where Rap, are winding resistance of stator winding, Lab, care self-inductance of each phase of the stator 
winding, and Lm is mutual inductance in any two phases. The back electromotive force is represented in (8), 
as a function of the rotor position @ at flux A, and speed w,.. It is worth noting that f (0) has the same trapezoidal 
shape as the counter-electromotive force per phase [20], where x represents the three phases (a, b, or c). 


ex = f O)Ax@r (8) 


The BLDC motor's output electromagnetic torque Temc is determined from the output power P, and the 
mechanical rotor speed is œw, [13], provided in (9): 


=, Ps = È exix (9) 


T „== 
emc wr wr 


The BLDC motor in our proposed system is connected to a centrifugal pump via its shaft [18]; the 
choice of this type of pump is related to its ability to operate at its rated speed and power for long periods of 
time even in the presence of partial shading, and thus its ability to pump a large volume of water while 
delivering a relatively high overall system efficiency. Water pumping necessitates a torque Tp proportional to 
the square of the motor speed w,. The characteristic torque looks (10) [19], [21]: 


Tp = Kp,” (10) 


3. MAXIMUM POWER POINT TRACKING CONTROL TECHNIQUES 
3.1. Particle swarm optimization 

In the search space of a PSO system, flies a swarm of individuals (particles) that represent a candidate 
solution to the optimization problem [22]. Each particle is fully characterized in terms of position z; and 
velocity @; as shown in (11) and (12). Initially, the particles are randomly initialized in the search space of the 
problem. Then, at each iteration, each particle moves and uses two final values to update their positions and 
velocities [23], [24]. The flowchart of the MPPT-based PSO is given in Figure 8(a). 
ghee het (11) 


i l 
it = wo + min (Ppest — Zi) + Marr (Gpest — 2°) a 


Where, w is called the inertia weight, m, and m, considered as acceleration coefficients, r4 and rz as random 
variables uniformly distributed within [0, 1], k is the iteration number, Pþest is as the personal best position of 
particle, and Gpest is called the best position of the particles for entire swarm population [23]. 

The reference voltage and output power of the PV system are determined sequentially by the position 
and velocity of the particle swarm when using the PSO approach in the MPPT domain [22]. As a result, (11) 
may be changed as (13): 

Vrepi = Vrega + Poa (13) 

Indeed, after each cycle, the particle with the greater power remains fixed in its position, while the other 
particles move in accordance with the position xi. Finally, after the required global maximum has been obtained, the 
ideal voltage is saved and returned as the output signal to be tracked PWM. The PSO optimizer runs as long as the 
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load requires power [25]-[27]. Clearly, changes in air conditions have an impact on power production and voltage. 
As a result, as (14) and (15) is used to validate the detection of fluctuations in power and voltage [27]: 


es > AP sot (14) 
\Ver.i — Vrep,a—1)| < Aset (15) 


Where AP,et and AV,,,; are the power and voltage thresholds, respectively. These ones are predetermined by 
the user (i), (i-1): the current and the previously measured values, respectively. The parameters of the PSO 
algorithm employed in our study are w=0.1, m,=m,=1.2, ray =raz=0.5, size of population is 5, and dimension 
number is 1. 


3.2. Grey wolf optimization 

GWO is a novel meta-heuristic algorithm inspired by grey wolves. The wolves in the pack are 
classified into four types based on their physical shape and strength: alpha (a), beta (B), delta (ô), and 
omega () [5], [28]. The GWO utilized (16)-(19) [29] and the flowchart of the proposed GWO-based MPPT 
algorithm is shown in Figure 8(b). 


D=|C.X,G) - X,Q)| (16) 
XG +1) =X,() — A.D (17) 
A=2ab,-a (18) 
C =2.b, (19) 


Where j is the current iteration; D, A, and C are the coefficient vectors; Xp is the position of the prey; X is the 
position of grey wolf; a is a variable which linearly decreased from 2 to O over course of iterations that 
resembles approaching the prey, and b,, b, are a random variable whose value is [0-1]. 
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Figure 8. Flowchart for (a) PSO MPPT algorithm and (b) GWO MPPT algorithm 
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When using the GWO approach to track maximum power, the reference voltage V„ep is specified as a 
gray wolf in this paper as shown in (20). Firstly, it is necessary to initialize the population (number of wolves) 
in our case it is considered three (the number of PV module) by setting minimum and maximum limits, then 
we calculate the PV power that determines the fitness value of the wolves, by determining two values as first 
and second maximum power. The search ends when the maximum PV power is obtained (the prey stops 
moving) after updating the positions. Each change of irradiation allows to reinitialize the search by validating 
as (21) [28], [30]. 


Vie G + 1) = Veep Q) — A. D (20) 
een (21) 


Where P, is the current global MPP of the panel. In the case of PSC, P,, will be much lower than Pm and AP 
is set at 0.05 to have better performance. 


3.3. Sliding mode controller design 

The SMC is a type of controller applied to non-linear systems, known for its robustness and ease of 
implementation compared to other non-linear control techniques [31]. A typical SMC includes two operating modes 
(bisectional). First, the sliding surface must be chosen in such a way that the system dynamics are stable, and then a 
control must be designed to make the switching surface appealing to the system state [27], [32], [33] as (22): 


y = Vv 
(22) 
A = Viet 
The tracking error € is (23): 
E =Y — Yréf (23) 
First, the sliding surface defined as (24): 
= d r-1 
=(w+<)e (24) 
With 'P is the positive constant, and r is the relative degree E. The derivative of output y is (25): 
; . 1 1 
y= Vpv = c, lev = rm (25) 
Using (3) and (6), the second derivative of y can be represented as (26): 
Per 1; 1 
ya Vov — c, lev = tc, Vw = (1 m d)V,) (26) 


As can be seen, the second derivation of y results in the appearance of the law control, which is the 
duty cycle d of the DC/DC boost converter. As a result, the relative degree r would be 2. Therefore, the sliding 
surface displayed in (24) will be (27): 


S=(¥+ Je (27) 


VA 
dt 
Its time derivative is expressed as (28): 
S=WE +Ë (28) 


Where the second time derivative is expressed as (29): 


Ë= y-y réf— Ve a Vref (29) 
By substituting (5) for (29): 


1: 1 + 
ë = Spv — 7 Vov — A - Ds) — Ürer (30) 
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Assuming Vier and Vier set to zero, the time derivative of the sliding mode surface is (31): 


ipv — —(Vpv — (1 — dK) (31) 


S= Mov + PY Tc 


The stability of the system can be obtained via the analysis of the Lyapunov theory defined positively as (32): 
V =<S? (32) 


It is only necessary that the derivative with respect to time of the Lyapunov function is negative while the 
coefficient k must be positive (33): 


V = SS = —ksign(S) (33) 


The convergence to zero in finite time of the sliding surface assures that the system will be asymptotically 
stable, the first-time derivative of S should be (34): 


1: 


; , ; 1 
$ = —ksign(S) = Mov + = low — iz V —(1-d)h,) (34) 
The duty cycle control d can be found as (35): 


C. r " ; V: v 
d= [ksign(S) + Vp» + Alyy | — EP +1 (35) 


4. SIMULATION RESULTS AND DISCUSSION 

The setup under consideration in this study consists of three PV modules linked in series, each with 
20 cells as shown in Figure 9, the selected module is AMETEK. MATLAB/Simulink was used to analyze the 
performance of the proposed system discussed in section 2. Table 1 shows the SPV specification and Table 2 
shows the boost converter parameter that was utilized to power the BLDC motor-water pump (TETRA 
142TR12 series) are defined in Table 3. Two irradiation profiles have been considered in this study: profile 1 
is the three PV modules receive initially between [0,1 s] an irradiation of 800 W/m?, 700 W/m?, and 
1,000 W/m? at 25 °C and profile 2 is between [1 s, 2 s] they receive 560 W/m?, 490 W/m?, and 700 W/m? at 
25 °C. It should be noted that these shifts are sudden, which is not realistic in practice. 


Model of One PV Module 


PV cells 
1-20 


PV cells 
21-40 Variable 
DC Source 


PV cells 
41-60 


Figure 9. Model of partially shaded PV panels 
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Table 2. The three-phase inverter’s switched states 


Hall signals Electronic commutation 

Intervals 
hı h, h} Sı S2 S3 S4 Ss s 
NA 0o o0 0 0 0 0 0 0 0 
0-60 0 1 1 0 0 oO 1 1 0 
60-120 1 1 0 0 1 1 0 0 0 
120-180 1 0 1 0 1 0 0 1 0 
180-240 1 0 1 1 0 0 0 0 1 
240-300 1 1 0 1 0 0 1 0 0 
300-360 0 1 1 0 0 1 0 0 1 


Table 3 .Parameters of BLDC motor 


Parameters Values 
Stator phase resistance 23 Q 
Stator phase inductance 0.05e-3 H 
Rated power 200 W 
Rated speed 600 rpm 
Voltage constant 11 V/krpm 
Flux linkage 0.052521 V/s 
Inertia 0.0008 kg.m? 
Viscose dumping 0.001 N.m.s 
Pole pairs 1 
Torque constant 0.10504 N.m/A 


4.1. Performance of solar PV array 

The two irradiance profiles considered in this paper are shown in Figure 10(a). A strong ripple in MPP 
voltage and power can be observed when PSO and GWO MPPT techniques are implemented, as shown in 
Figures 10(b) and (c), respectively. The PSO MPPT technique expects a voltage of 32.95 V, while that achieved by 
GWO is 32.66 V. The maximum power seen in the two MPPT methods proposed in this work is 186.3 W (PSO and 
GWO). Because the system under consideration begins with a radiation of 700 W/m”, there are high transients 
induced by the GWO method at all levels of irradiation in the profile 1, resulting in a significant loss of power. When 
the shade pattern changes in profile 2, a new P-V curve is depicted in Figure 10(b), and both algorithms search the 
P-V curve for a new MPP. The combined maximum power of the two systems shown is 104 W. The PSO method 
allows for a 0.3 s stabilization time to attain the MPP in the transient state during [0-1 s]. The graphs of voltage and 
PV power are zoomed in on the 0.2-1 s intervals to highlight the transient and steady state ripples, respectively. 
During the interval [1,2], both methods exhibit the same performance and voltage and maximum power levels. As 
the irradiance changes, the SPV array's power tracking improves dramatically. Under shifting climatic 
circumstances, the MPPT approaches clearly optimize the SPV array power. 


Cl 
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Figure 10. Performance of SPV array under PSC: (a) irradiance profile, (b) PV generator voltage, and (c) PV 
generator power 
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4.2. Performance of DC-DC boost converter 

The performance of the boost converter is shown in Figures 11(a) and (b). They are associated with the 
duty cycle in Figures 11(c) and (d). The continuous conduction mode of operation is maintained even when the 
irradiance varies. Boost's conductance current J; is 5.8 A when PSO MPPT is used, and it increases from 4.72 A to 
5.76 A when GWO is used for [0-1 s]. The J, conductance current is 4.22 A for the second profile between [1-2 s], 
according to both MPPT algorithms investigated. The DC link voltage is represented by the boost output voltage V; , 
varies from a low of 23 V during [1-2 s] to a high of 32 V during profile 1 [0-1 s]. When compared to the GWO 
approach, the PSO algorithm greatly reduces the ripple content of the inductor current and link voltage. 
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Figure 11. Boost performance when implementing the PSO and GWO MPPT techniques: (a) link 
voltage V (s), (b) inductor current IL, (c) duty cycle for PSO MPPT, and (d) duty cycle for GWO MPPT 


4.3. Performance of BLDC motor-pump 

Figures 12(a)-(d) (in Appendix) for GWO MPPT and Figures 13(a)-(d) (in Appendix) for PSO MPPT 
show the BLDC motor parameters. In the first profile, the BLDC motor reaches its maximum speed of 
522 rpm in 0.4 s using PSO MPPT, while the GWO algorithm require 0.5 s. Due to electronic commutation, 
the BLDC motor creates electromagnetic torque with pulsations. The PSO MPPT method generates the least 
torque ripple. The soft-switching technique regulates the motor stator current via the DC-DC boost, and it 
achieves the rated intermediate circuit current of 5.78 A in 0.08 s. 

The PSC slows down the pace in the second profile. The waveform in Figures 11(c) and 13(c) shows 
that the BLDC motor achieves 399 rpm using both approaches given. The motor may also run at 50% of the 
irradiance, or 500 W/m?. As illustrated in Figures 12(d) and 13, the speed and electromagnetic torque may 
smoothly follow irradiance fluctuations. Due to the electrical switching of the VSI, pulses can be seen in the 
electromagnetic torque of the BLDC motor Figure 12(d). 


5. CONCLUSION 

This study describes how to use SMC theory on a DC-DC boost converter to power a water pump 
with a BLDC motor under various climatic circumstances. The flexible control provided by PSO and GWO 
MPPT has enabled the intelligent and efficient operation of the solar water pumping system via the boost to 
satisfy the initial objectives of MPPT extraction, smooth speed variation, and soft starting of a BLDC motor. 
The designed controller was able to predict the duty cycle for the boost converter to extract the maximum 
available PV power and ensure the soft start of the BLDC motor. The transient and steady state performance 
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of the PSO MPPT method compared with GWO reveals that the power tracking efficiency of the PSO MPPT 
algorithm is the highest. 
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Figure 12. Performance of BLDC motor-pump under variable solar irradiation with GWO algorithm: 
(a) back emf ea, (b) stator current I4, (c) stator speed w,., and (d) electromagnetic torque Te 
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Figure 13. Performance of BLDC motor-pump under variable solar irradiation with PSO algorithm: 
(a) back emf ea, (b) stator current J, 
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Figure 13. Performance of BLDC motor-pump under variable solar irradiation with PSO algorithm: 
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